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Resumen 


Se desarrolló un procedimiento para dimensionar zanjas de infiltración dependien- 
do de un balance hídrico en la ladera, el cual está focalizado en dos elementos 
centrales: eventos extremos de precipitación y capacidad de infiltración constante. 
Debido a la disponibilidad de datos hidrológicos de la cuenca del río Purapel, 
este modelo de infiltración incorpora análisis regional de frecuencia, análisis de 
distribución de la lluvia, método de índice de tormenta y técnicas geoestadísticas. 
Los resultados permiten obtener una extensión de las relaciones intensidad-duración 
frecuencia para zonas sin registro continuo de la lluvia y el dimensionamiento para 
la estructura de conservación de aguas y suelo. El modelo fue aplicado a la cuenca 
del río Purapel, en Chile central, para demostrar la facilidad del procedimiento 
sugerido. 
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Introducción 


El marcado impacto de la erosión hídrica en 
la disminución de la productividad agrícola 
y forestal en la zona central de Chile hace 
necesario desarrollar planteamientos y diseños 
de obras de conservación de aguas y suelos. 
Así, cabe destacar que la agricultura de 
escorrentía y la agrosilvicultura son medidas 
de aprovechamiento hídrico para un cultivo, 
con favorables resultados de establecimiento, 
crecimiento y desarrollo de plantas en zonas 
semiáridas. Este trabajo estudia las zanjas de 
infiltración, para aumentar la recarga hídrica 
del suelo (Boers y Ben-Asher, 1982; Lovenstein 
et al., 1991). Se han propuesto algunos modelos 
interesantes (Falco et al., 1997; Martínez de 
Azagra, 2000; Akan, 2002), que suponen que 
toda el agua de escorrentía generada por el área 
de impluvio se infiltra en el área de recepción, 
situación que dista mucho de lo real. 


En muchos de estos modelos se dispone de 
escasa información climatológica. Además hay 
una gran variabilidad espacial y temporal de los 
procesos hidrológicos. Ello requiere incorporar 
análisis regionales y técnicas geoestadístiscas. 
(1998) propusieron 
unas curvas IDF regionales que permiten la 


Koutsoyiannis et al. 


incorporación de estaciones sin registros de 
pluviometría (con distintas duraciones de 
lluvia) para la construcción de la IDF en sitios 
no aforados. La determinación de un índice 
de tormenta derivado de las relaciones IDF 
permite establecer patrones de las intensidades 
de precipitación para diferentes duraciones de 
lluvia (Di Baldassarre et al., 2005). La variabilidad 
de este índice de tormenta se caracteriza 
mediante métodos estadísticos para predecir la 
distribución espacial de los valores de la variable 
intensidad de precipitación, expresado como un 
cociente o índice de tormenta, 1”, / *,, (Cheng ef 
al., 2003). También estos modelos requieren de 
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Cuadro 1. Modelos intensidad-duración-frecuencia (IDF). 
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Modelo IDF T,= Número de parámetros 
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una determinación de las propiedades físicas 
del suelo, en especial la capacidad final de 
infiltración del suelo, f.. 

La opción más sencilla es trabajar con 
eventos de intensidad de lluvia y capacidad 
de infiltración constante, aun cuando se 
podría subestimar el calado (1) almacenado 
y sobreestimar la superficie de la infiltración 
(Akan, 2002). 

Este trabajo tiene como objetivo desarrollar 
un modelo de zanja de infiltración en función 
del balance hídrico de una ladera usando un 
sistema de información geográfico. 


Metodología 
Diseño hidrológico de zanjas de infiltración 


El diseño de un sistema de aprovechamiento 
de agua (figura 1) se define como un sistema 
“producción-recolección” en una ladera que 
se subdivide en N unidades. Las hipótesis del 
modelo de zanjas de infiltración son: 


e La precipitación es uniforme en cada 
unidad analizada. 

+ Se considera un suelo desnudo. 

+ No hay evaporación desde el suelo durante 
el chubasco. 

+ Las unidades son independientes entre sí 

e El tamaño de cada unidad < 100 tm?. 

+ No hay circulación de agua entre las 
distintas unidades. 


e La sección transversal y el calado son 
uniforme a través de la zanja de infiltración. 

+ Ante eventos prolongados de lluvia y el 
parámetro a y f sea grande, la velocidad 
de infiltración es: 


zona de recepción: o, (1) =£+(f-f.)Je7* de 
simplifica a v, (+) =f. 

zona de impluvio: w,(t)=3.+(2)-3,Je Y, se 
simplifica a vw, (+) =2. 


siendo v(t) = velocidad de infiltración de 
la zona de recepción; w(t) = velocidad de 
infiltración en el área de impluvio; f, y g, = 
capacidad de infiltración inicial del área de 
recepción e impluvio, y f. y 2. = 
de infiltración final del área de recepción e 


capacidad 


impluvio, respectivamente. Por lo general f = g. 
(salvo enmiendas muy importantes en el área 
de recepción o aplicaciones impermeabilizante 
en el área de impluvio). 

Este diseño plantea que el volumen aportado 
(V,) sea igual al volumen recepcionado (V). La 
determinación del volumen de aportación se 
basa en el cálculo de la escorrentía superficial, 
explicado por la ecuación de continuidad: 

q=P-P,-F (1) 

Donde q es la escorrentía superficial [L]); 
P, la precipitación total [L]; P,, la abstracción 
inicial [L]; F, la infiltración acumulada [L]. 
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Figura 1. Esquema de zanjas de infiltración en ladera. 


La precipitación P, en un intervalo At, puede 
expresarse usando las relaciones intensidad- 
duración-frecuencia (Chen, 1983): 


: (2) 
t+d 


siendo 1, la intensidad media de precipitación 
[LT"]; t, la duración del evento [T]; T, el periodo 
de retorno [T]; a, b, c y d, parámetros del modelo 
IDF. Derivando la ecuación (4), se obtiene la 
intensidad de lluvia instantánea, i, [LT"]: 


(3) 


Si se considera una superficie suficien- 
temente pequeña, se puede infiltrar toda el 
agua. Transcurrido un tiempo f = É, (tiempo 
de encharcamiento), el suelo alcanzará su 
estado de saturación y su superficie comenzará 
a encharcarse, llenándose las depresiones. 
A partir de este instante, el sistema produce 
escorrentía, que termina cuando la intensidad 


de lluvia instantánea (i,) se iguala a la capacidad 
de infiltración final (f.), de modo que: 


(4) 


El valor de t correspondiente a q (ecuación 
(1)) puede ser obtenido usando el método de 
Newton-Raphson. 

En zonas geográficas con disponibilidad de 
datos pluviométricos, la intensidad media de la 
lluvia puede ser fácilmente obtenida utilizan- 
do la curva local de intensidad-duración- 
frecuencia (Yu y Chen, 1997). Sin embargo, esto 
no ocurre en zonas de la región central de Chile, 
por lo que hay que recurrir al análisis regional 
(Hosking y Wallis, 1997), análisis de distribu- 
ción de precipitaciones máximas (Raynal, 
2005), índice de tormenta (Di Baldassarte et al., 
2005) y técnicas geoestadístiscas (Cheng et al., 
2003), como el krigeado ordinario. La lluvia 
de diseño se determina a partir del uso de un 
índice de tormenta, k: 

TÚ =kxll 


diaria 


(5) 


Tecnología y 


Ciencias del Agua, vol. IL, núm. 2, abril-junio de 2012 


Ciencias del Agua, vol. IL, núm. 2, abril-junio de 2012 


/ 


(vy80]0492], 


Flores-Villanelo, Diseño de zanjas de infiltración en zonas no aforadas usando SIG 


donde I', es la intensidad de precipitación a la 
hora í (i=1, 2, 6, 12 y 24 horas) e I” 


diaria 


es la 
intensidad de precipitación diaria. 

La substracción inicial, P, [L], se determina 
por Ferrér et al. (1995): 


100 ) (6) 


2, =508 [71 
CN 


siendo CN el número de curva para condicio- 
nes de humedad antecedente del suelo tipo III. 
La infiltración acumulada es: 
F=f.-t,/60 (7) 
Donde t,=t-+t,, es el periodo de infiltración 
[T] y t,, =(P,/i,)*60, el tiempo que corresponde 
a la substracción inicial [T]. 
Para el dimensionamiento de una zanja de 
infiltración rectangular en ladera se consideran 
los siguientes volúmenes de agua: 


Volumen de agua de impluvio: 


V, =9-S¡=q-1-D (8) 
Volumen geométrico de la zanja: 
EG =A75, 470,1 (9) 


Volumen de precipitación al interior de la zanja: 


V, = (1, t)S, =P-b,:1 (10) 
Volumen de infiltración: 
Y =(0,:t):S, =(£:t):b,:1 (11) 


Siendo D el distanciamiento entre zanjas 
[L]; h, la profundidad de la zanja [L]; S, el área 
de recepción dentro de la zanja [L?]; b, el ancho 
de la zanja [L]; !, el largo de la zanja [L]; v, la 
velocidad de infiltración de agua dentro de la 
zanja (f.) [LT"]. 

Utilizando la ecuación de continuidad apor- 
tación-recepción (V, = V/): 


q:S =h¿S,+ f,-t:S,-P-S, (12) 


sujeto a h, >t(£. = in) 


se deduce el espaciamiento entre zanjas, D [L], 
y con ella la del intervalo vertical entre zanjas 
de infiltración, í, [L], para aplicaciones de un 
SIG, con el modelo de digital de elevaciones de 
una cuenca, por lo que se tiene que: 


D=b,[h, +t(f. -i,)](gcos(o))” (13) 


Área de estudio 


El área de estudio (figura 2) comprende 
la cuenca del río Purapel, de clima tipo 
mediterráneo, con seis meses secos y seis 
meses húmedos, ubicada en la cordillera de 
la costa de la región del Maule, Chile. Cuenta 
con un área de tributación de 264.6 km? hasta 
la estación de aforo en Nirivilo (35% 34' LS; 72” 
05' LO), con un largo de cauce principal de 77.2 
km. La temperatura media anual es de 15.2 *C, 
con temperaturas máximas promedio de 28.5 
“C en enero y mínimas promedio de 5.4 *C en 
julio. Las precipitaciones medias anuales en la 
cuenca se estiman en 810 mm, siendo julio el 
mes más lluvioso, con 194.2 mm, y enero el mes 
más seco, con 4.4 mm (Pizarro et al., 2006). 

Las precipitaciones en la zona central de 
Chile ocurren en un 90% durante el periodo 
mayo-septiembre, asociadas con los frentes 
fríos que ocasionalmente alcanzan esta región 
cuando el anticiclón subtropical del pacífico 
suboriental es anormalmente débil (Garreaud y 
Rutllant, 2006). Así, esta zona se caracteriza por 
una extrema variabilidad de las precipitaciones 
entre un año y otro, en parte relacionada con la 
Oscilación del Sur de El Niño (ENSO), mientras 
que el comportamiento interdecadal se asocia 
con la oscilación decadal del pacífico (PDO) 
(Montecinos y Aceituno, 2003). La clase textural 
predominante del suelo es franco areno-arci- 
llosa y tiene un comportamiento hidrológico 
correspondiente al grupo C de acuerdo con la 
clasificación del método del número de curva. 
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La cuenca presenta una cobertura forestal 
sobre el 70%, entre plantaciones de Pinus 
radiata (D. Don), actualmente en explotación, y 
bosque nativo. En los terrenos más planos de la 
cuenca se presentan cultivos característicos de 
la agricultura de secano. 


Determinación del diseño hidrológico de zanjas 
usando ArcGis 


La estimación de la lluvia de diseño se deter- 
minó mediante la ecuación (5). Se analizó 
la presencia de datos incorrectos y valores 
anómalos, independencia y estacionariedad, y 
tendencia monótona de las series (1965-2007), 
a través del test de Grubbs y Beck, Wald- 
Wolfowitz y Mann-Kendall, respectivamente. 
El análisis regional se llevó a cabo con los 
programas XTEST, XSIM y XFIT, descritos 
ampliamente en García (2000), como una 
estimación de observaciones puntuales de 
las precipitaciones máximas de la región 
central de Chile. Se hizo necesario interpolar 
la información suministrada por los 52 
pluviómetros de la misma, mediante krigeado 
ordinario. 

El proceso geoestadístico para la variable 
“índice de tormenta, k” ajustó un modelo 
exponencial a una serie de 52 puntos (o 
estaciones pluviométricas) según Isaaks y 
Srivastava (1989): 


ha 


y(h1¿)=Cj +C; leal :) cuando |, >0 (14) 


y((1,)=0 cuando lh,]=0 (15) 


Donde C, es la varianza en el origen; 2, 
el alcance; h,, el vector distancia, y C, + C,, la 
meseta. 

Los parámetros del variograma se consi- 
deraron funciones también del intervalo 
de recurrencia y de la duración de la lluvia 


de diseño. Por lo tanto, se determinó un 


variograma para cada diseño de lluvia. 
Esas relaciones son bien consistentes con las 
características aleatorias de la lluvia de diseño, 
considerando que la distribución espacial de la 
lluvia tiene un comportamiento (influenciado 
por el anticiclón del Pacífico) como un campo 
aleatorio y su estructura de variación espacial 
varía en función de la duración y el periodo de 
retorno. 

Los valores de intensidad de precipitación 
interpolados a la cuenca del río Purapel son 
ajustados a un modelo IDF mediante el méto- 
do de mínimos cuadrados y utilizando el 
algoritmo de optimización de Rosenbrock, a 
través del método de gradiente avanzado de 
quasi-Newton. Se probaron varios modelos. 

Este modelo estima de forma indirecta la 
velocidad de infiltración del agua en el suelo 
(figura 3), en función de pendientes y textura 
del suelo, para lo cual se incorporó la base 
de datos de la caracterización de suelos del 
estudio agroecológico (1:50.000) del Centro de 
Información de Recursos Naturales (Ciren)- 
Chile y las tasas de infiltración constante 
publicadas por el Departamento de Suelos de 
la Universidad de Carolina del Norte, Estados 
Unidos. 

La aplicación SIG de este modelo se lleva a 
cabo mediante un código Matlab y un modelo 
de elevación digital de la cuenca, obtenido de 
la información digitalizada de una ortofoto 
del Centro de Información de Recursos 
Naturales 1:20.000 (Ciren, Chile) y cartas 
digitales del Instituto Geográfico Militar 
1:50.000 (IGM, Chile), en donde se utilizó 
ArcGis de ESRÍ para el ordenamiento de la 
información disponible y la visualización 
gráfica de ésta. 


Resultados y discusión 


La exploración de los datos de precipitación 
extrema registró un 1% de valores anómalos 
del total de datos de las series, donde el 21% 
de los mismos se registró en la serie de una 
hora, por lo cual se eliminaron del análisis 
regional de precipitaciones. Este estudio 
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Región de Coquimbo 


Región de Valparaíso 
Región Metropolitana 
(Santiago de Chile) 


Región Lib, Bernardo 
O'Higgins 


Región del Maule 


Región del Bío-Bío 


Región de la Araucanía 


Figura 2. Ubicación de la cuenca del río Purapel y de las 52 estaciones pluviométricas de la región central de Chile. 


consideró un nivel de significancia del 1%, 
para lo cual la totalidad de las series de 
precipitaciones máximas cumplen el requisito 
de independencia y estacionariedad. En 
cuanto al análisis exploratorio de la serie 
Nirivilo, la prueba de Grubbs y Beck sí 
detectó un dato anómalo (en un nivel de 
significancia del 10%), en la lluvia del día 26 
de abril del año 1957, que se eliminó. Mientras 
que el estadístico de  Wald-Wolfowitz 
comprobó que los datos de precipitaciones 
son independientes y estacionarios, y la 
prueba no-paramétrica de Mann-Kendall no 
detectó tendencias monótonas en la serie. 


Los estimadores de los parámetros para la 
distribución GVE, por medio del método 
de momentos de probabilidad ponderada 
propuesto por Raynal (2005), fueron muy 
cercanos en valor numérico a la metodología 
de Hosking y Wallis tanto en la obtención 
de parámetros de la función de distribución 
GVE como en la evaluación de los valores de 
diseño, sin embargo, el insesgamiento de los 
primeros es menor. 

Los modelos matemáticos ajustados 
presentan excelentes resultados (mayoritaria- 
mente R?> 90%). Para la cuenca del río Purapel, 


el mejor fue el modelo 1 (Témez, 1978): 
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Textura 

0-3% 
Arena 20 
Franco arenoso 1.8 
Franco arenoso fino 1.3 
Franco y franco limoso 0.8 
Franco arcillo arenoso y franco arcillo limoso 0.5 
Arcilla, arcilla arenosa y arcilla limosa 0.25 


Tipos de suelo presentes en la cuenca 


Asociación Alto Colorado 


Pendiente 
> 3-9% > 9% 

1.8 1.3 
1.3 1.0 
1.0 0.8 
0.5 0.4 
0.4 0.3 

0.15 0.08 


Asociación Constitución 


¡EE Asociación Treguaco 
MA Carampangue 
EA] Misceláneos 

MEA] Ninhue 


Figura 3. Tasas de infiltración aproximadas (cm /h) para varias texturas y pendientes del suelo para la cuenca 


del río Purapel (fuente: http: / /www.soil.ncsu.edu). 


r  19.37-T02109 


NES 5 (16) 
E 03501 ,0,5677 


Por otra parte, distintas relaciones de f. en 
función de la duración de la lluvia (t); periodo 
de retorno (7); escorrentía superficial (q); y 
distanciamiento entre zanjas de infiltración se 
muestran en las figuras 4, 5, 6 y 7. En la figura 
4, el tiempo que demora en igualar las tasas 
de intensidad de precipitación instantánea (i,) 
y la capacidad de infiltración (f.) disminuye 
al aumentar fc, y demuestra que el modelo 
IDF 1 se comporta mejor que los modelos 4 
y 5, propuestos por Yarnell (1935) y Ponce 
(1989), respectivamente. La figura 5 presenta 
una familia de curvas para cada periodo 
de retorno (5, 10 y 20 años) para la relación 
t vs. Le Al aumentar el periodo de retorno, 
el tiempo de equilibrio (i = f) aumenta. 


Por otra parte, la escorrentía superficial 
disminuye con f. (figura 6); contrariamente, 
el distanciamiento entre zanjas aumenta 
con f. (figura 7). Para periodos de retornos 
altos, D comienza a disminuir debido a que 
los volúmenes de escorrentía a infiltrar son 
mayores, para lo cual se deben realizar más 
obras de infiltración por hectárea. 

El modelo DEM de la cuenca del río 
Purapel (figura 8) es una red de 247 x 266 
celdas, incluidos pendientes, cauces hídricos y 
retículas de banco. El tamaño de la celda fue 
definido como 100 * 100 m (1 ha). Al convertir 
las capas de información ASCH a formato 
matricial, provenientes de Matlab, el modelo 
representa en ArcGIS el distanciamiento entre 
zanjas de infiltración para cada celda dentro de 
la cuenca (figura 9). 

El distanciamiento entre zanjas está 
relacionado con una serie de variables y 
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Figura 4. Duración de la lluvia cuando ¡,= f. y la capacidad Figura 6. Escorrentía superficial en función de la capacidad 
de infiltración (f.) del suelo, para T=10, usando tres de infiltración del suelo (T = 10; 8 = 10). 
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El Figura 5. Duración de la lluvia cuando i,=f. y la capacidad Figura 7. Distanciamiento entre zanjas en función de la 
E de infiltración (f.) del suelo, para distintos periodos de capacidad de infiltración del suelo (T = 10; 0 = 10; 
Ñe retorno (T=5, 10, 20). b=0.5;h=0.3). 
£ 
E 
S parámetros, los cuales el ingeniero de obra obra, h, provoca un mayor D (figura 10). Es 
ed puede modificar según los objetivos del posible considerar un factor de rebosamiento, 
E proyecto de diseño. Aumentar la altura de la Je lo cual permite manejar D, cuando 
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El 175-230 
E 231-285 
[1 286-340 
A] 341-395 
A 396-450 
MM 451-505 
A 506-560 
MM 561-615 
E 616-670 
MM 671-725 


Altitud 
msnm 


Pendientes (%) 


El 0-3.2 
El 3.3-8.1 
El 8.2-12.3 
HA] 12.4-16.0 
Má 16.1-19.6 
MA 19.7-23.2 
MÁ 23.3-27.0 
MM 27.1-31.3 
MM 31.4-36.6 
MM 36.7-48.3 


Figura 8. Modelo digital de elevaciones de la cuenca del río Purapel (izquierda) y mapa de pendientes de la cuenca del río 


Purapel (derecha). 


(a) Escorrentía 
superficial (mm) 


OJo0<0<3f=3.81 
Mi 3<0<9f =7.62 
Mi  9>9f=7:52 


(b) Capacidad de 
infiltración (mm /h) 
textura: SCLm 


(c) Distanciamiento 
zanjas (m) 


[El 3.5-4.6 
HA 4.7-5.6 
Ll 5.7-6.7 
EX 6.8-7.7 
MM 7.s-8.7 
MM s.s-9.s 


Figura 9. (a) Mapa de escorrentía superficial (izquierda); (b) capacidad de infiltración del suelo f. para una textura franco areno- 


arcilloso (centro), y (c) distanciamiento entre zanjas de infiltración (derecha), para T =10; b =0.5;h=0.3, aplicado a la cuenca 


del río Purapel. 


los volúmenes de impluvio exceden las 
necesidades del cultivo. Al aumentar el factor 
de rebosamiento, se incrementa D (figura 11). 
Otras relaciones pueden encontrarse (figuras 12 
y 13) al relacionar distintos valores del número 
de curva y la pendiente, y el distanciamiento 
entre zanjas, manipulando la configuración 
geométrica de la zanja (en este caso, h). 

Para comprender las relaciones entre D 
y CN II (figura 12), hay que considerar que 
el suelo está desnudo; de lo contrario debe 
incorporar nuevos procesos en el modelo, 


como la intercepción y, posiblemente, la 
evaporación. 

Uno de los aspectos más discutibles del 
presente modelo es la variabilidad temporal 
y espacial de las variables incluidas, como 
la velocidad de infiltración constante en 
el impluvio y recepción, la intensidad de 
precipitación constante, la superficie anegada 
o el volumen de agua almacenado en las 
microdepresiones, etcétera. Así, el valor de 
q y h disminuye cuando f aumenta. Este 
comportamiento puede ser bien entendido 
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Figura 10. Distanciamiento (D) en función de la altura de la 
zanja de infiltración (1) para distintos periodos de retorno 
(T =2, 5, 10, 20; 0 =10; b =0.5). 
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Figura 12. Distanciamiento (D) en función del número de 
curva (condición hidrológica III) para distintas alturas de 
zanja de infiltración (h = 0.3, 0.4, 0.5, 0.6; 0 = 10; b =0.5). 
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Figura 11. Distanciamiento (D) en función de la altura 
de la zanja de infiltración (/1) para distintos factores de 
rebosamiento (f, = 45 / R(ninguno), 0.9 (10%), 70.7 (30%), 
90.5 (50%); 0 =10; b =0.5). 


cuando se revisan los resultados de Pruski ef 
al. (1991), quienes obtuvieron gráficas para q 
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Figura 13. Distanciamiento (D) e intervalo (i,) entre zanjas, 
en función de la pendiente del terreno (0) para distintas 
alturas de zanja de infiltración (h = 0.3, 0.4, 0.5, 0.6; 8 = 10; 
b=0.5). 


con la duración de la lluvia y diferentes valores 
de infiltración (f). Al aumentar f, la duración 
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de la lluvia que produce q disminuye, como 
también el valor total de la lluvia asociada 
con éste, puesto que como la duración de la 
lluvia aumenta, estas intensidades máximas 
promedio disminuyen, mientras que el total de 
lluvia se incrementa. 

Asimismo, al no validar los coeficientes del 
parámetro P, (ecuación 6), pueden limitar los 
resultados, ya que estos fueron desarrollados 
para Otras zonas geográficas, pudiéndose 
utilizar como herramienta de investigación de 
valores iniciales. Por ello, para el diseño óptimo 
de un sistema de captura y aprovechamiento 
hídrico, es necesario una determinación de alta 
precisión de los valores de f.. Para posteriores 
investigaciones se realizará una adaptación del 
número de curva a las nuevas fuentes de datos; 
en concreto, a la nueva cartografía de usos de 
suelo provenientes del proyecto Land Cover 
Corine (CLC-2000), ampliamente usado en 
Europa y Latinoamérica. 


Conclusiones 


El correcto diseño de una zanja infiltración para 
un periodo de retorno dado depende en gran 
medida de la precisión en la determinación de 
los volúmenes de escorrentía y la capacidad de 
infiltración (f.). 

Los resultados del análisis local y regional 
de las series anuales máximas de precipitación 
en la estación Nirivilo (cuenca del Purapel) y 
la región central de Chile, respectivamente, 
permiten obtener un valor confiable de I”, que 
reproduce además las propiedades estadísticas 
de las lluvias extremas observadas para la 
región de estudio. 

La utilización de un sistema de información 
geográfica y el krigeado ordinario resultaron 
ser una eficaz técnica de interpolación de datos 
muestrales y permite realizar una cartografía 
de índice de precipitación, k (1, /1%,,); asimismo, 
genera un estimador lineal e insesgado que 
hace mínima varianza del error de estimación. 
Esto permitió la construcción de las relaciones 
IDF para la cuenca del río Purapel, herramienta 
fundamental para el diseño de proyectos de 


conservación de aguas y suelos para zonas sin 
registros pluviométricos. 

Finalmente, el diseño de zanjas introducido 
en este trabajo emplea una base física mucho 
más extensa que otros diseños usados en 
la práctica; sin embargo, utiliza una mayor 
cantidad de información, lo que no siempre 
se dispone a cabalidad en la región central 
de Chile, y es que en la caracterización de 
los procesos de infiltración en una ladera con 
estructuras de 
posibilidades, desde trabajar directamente 


infiltración caben distintas 


con datos experimentales obtenidos sobre el 
terreno con infiltrómetros o con simuladores 
de lluvia, hasta utilizar modelos sencillos 
como el método SCS, o trabajar con tasas de 
infiltración constante, o bien, con modelos de 
infiltración en función del tiempo y volúmenes 
de agua precedentes. El diseño presentado 
aquí es uno de ellos y puede resultar una 
herramienta fundamental de planificación y 
diagnóstico para muchos proyectos que se 
pueden llevar a cabo en laderas “con procesos 
de desertificación” o de “establecimiento de 
cultivos agroforestales”. 
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Abstract 


FLORES-VILLANELO, J.P. Infiltration trench design for ungauged areas using SIG. Water 
Technology and Sciences (in Spanish). Vol. II, No. 2, April-June, 2012, pp. 27-39. 


A procedure was developed for sizing infiltration trenches based on the water balance of a 
hillslope. This is focused on two main elements—extreme rainfall events and constant 
infiltration capacity. Because of the shortage of hydrological data neat the site of the water 
resources project, this infiltration model uses regional frequency analysis, analysis of rainfall 
distribution, the storm index method and the application of geoestatistical techniques. The 
results allow for obtaining an extension of the intensity-duration-frequency relationships with 
non-recording rain gauges and the sizing ofwater and soil conservation structures. The model 
was applied to the Purapel River Basin located in central Chile to demonstrate the ease of the 
suggested procedure. 


Keywords: infiltration rate, intensity-duration-frequency, infiltration trench. 
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